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Abstract. We summarize our results on the Galactic merger rate of 
double neutron stars (DNS) in view of the recent discovery of PSR J0737— 3039. 

\ We also present previously unpublished results for the global probability 

J> ■ distribution of merger rate values that incorporate the presently known 

systematics from the radio pulsar luminosity function. The most likely 

■ value obtained from the global distribution is only ~ 15Myr _1 , but a 

!y-^ . re-analysis of the current pulsar sample and radio luminosities is needed 

' for a reliable assessment of the best fitting distribution. Finally, we use 

our theoretical understanding of DNS formation to calculate a possible 
upper limit on the DNS merger rate from current Type Ib/c supernova 

. rate estimates. 

6 ■ 
a: 

1. Introduction 

Soon after the discovery of the highly relativistic pulsar J0737— 3039 (Burgay et 
al. 2004) we applied our analysis method for pulsar populations and calculated 
the updated merger rate estimates for the current sample of Galactic close DNS 
(Kalogera et al. 2004). Our main conclusion was that this new, remarkably 
relativistic system is very significant for these estimates and leads to a rate 
increase by a factor of 5 — 7. This implies a correspondingly significant increase in 
DNS inspiral event rates for gravitational- wave (GW) interferometers like LIGO. 
In what follows, we summarize our recent results and present new results on: (i) 
how a global probability distribution of rate estimates can be calculated with the 
systematic uncertainties; (ii) the possible upper limits that could be imposed on 
the DNS merger rate if we adopt the theoretically expected relationship between 
DNS merger rates and Type Ib/c supernova (SN) rates. 



2. The revised Galactic DNS merger rate 

In Kim, Kalogera, & Lorimer (2003; hereafter KKL), we introduced a statistical 
method to calculate the probability density function (PDF) of the rate estimates 
for Galactic close DNS. After the discovery of PSR J0737— 3039, we derived a 
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combined P(TZ) considering the three observed DNS systems in the Galactic 
disk (for details see Appendix A of Kim et al. 2004). 

To calculate the merger rate of DNS systems in our Galaxy, we need to 
estimate: (i) the number iVtot of Galactic pulsars with pulse and orbital char- 
acteristics similar to those in the observed sample; (ii) the lifetime riif e of each 
observed system; (iii) an upward correction factor /b for pulsar beaming. 

We calculate 7V tot by modeling in detail the pulsar-survey selection effects for 
a number of pulsar population models described in KKL. The model assumptions 
for the pulsar luminosity function dominate the systematic uncertainties of our 
overall calculation. 

The lifetime of the system is defined by riif e = r sc j + r mrg , where r s d is a 
spin-down age of a recycled pulsar (Arzoumanian, Cordes, & Wasserman 1999) 
and r mrg is the remaining lifetime until the two neutron stars merge (Peters & 
Mathews 1963). We note that the lifetime of J0737— 3039 is estimated to be 
185 Myr, which is the shortest among the observed systems. 

The beaming correction factor is defined as the inverse of the fractional 
solid angle subtended by the pulsar beam. Its calculation requires detailed 
geometrical information on the beam. Following Kalogera et al. (2001), we adopt 
/ b = 5.72 for PSR B1913+16 (Hulse & Taylor 1975) and 6.45 for PSR B1534+12 
(Wolszczan 1991). Without good knowledge of the geometry of J0737— 3039A, 
we adopt the average value of the other two systems (~ 6.1). 

In Figure 1, we show P(7Z) for our chosen reference model that allows for 
a low minimum pulsar luminosity (Model 6 in KKL). The most likely value of 
1Z turns out to be 83 Myr -1 , larger by a factor of ~ 6.4 than the rate estimated 
before the discovery of J0737— 3039. We find the same increase factor for all 
pulsar population models examined. This revised merger rate implies an increase 
in the detection rate of DNS inspirals for ground-based GW interferometers such 
as LIGO (Abramovici et al. 1992). Using the standard extrapolation of our 
reference model out to extragalactic distances (see Kalogera et al. 2001), we 
find that the most probable event rates are 1 per 29yrs and 1 per 2 days, for 
initial and advanced LIGO, respectively. At the 95% confidence interval, the 
most optimistic predictions for the reference model are 1 event per 8yrs and 2 
events per day for initial and advanced LIGO, respectively. For more details see 
Kalogera et al. (2004). 

The revised DNS merger rate is dominated by PSR J0737-3039. Therefore, 
if the estimated lifetime of this system is revised in the future, it will directly 
affect our rate estimation. Lorimer et al. (in this volume) calculated the spin- 
down age of the system with various spin-down models and suggested an age 
in the range 30-70 Myr, which is shorter than the value we adopted for our 
calculation (r s d =100 Myr). The edges of this range give us rate estimates of 
TZ ~ 90 - 115 Myr -1 . 

The beaming correction for J0737— 3039 is also important for our rate es- 
timation. Recently, Jenet & Ransom (2004) suggested a geometrical model for 
this newly discovered system. According to their model, the predicted beaming 
correction factor is >6 assuming a two-sided beam (Jenet 2004, private com- 
munication). If confirmed, this could lead to a further dramatic increase in the 
merger rate estimates. 
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Figure 1. The PDF of DNS merger rate P(TZ) is shown on a log scale. 
The thick solid line is the total Galactic rate estimate overlapped with 
results for individual observed systems (dashed lines). Dotted lines 
indicate confidence intervals for the rate estimates. The same results 
are shown on a linear scale in the small inset. All results shown are for 
our reference model. 
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3. Global probability distribution of the rate estimates 

In KKL, we showed that estimated Galactic DNS merger rates are strongly- 
dependent on the assumed luminosity distribution function for pulsars. So far, 
we have reported results for each set of population model assumptions. Here we 
describe how we can incorporate the systematic uncertainties from these models 
and calculate, P g {1Z), a global PDF of rate estimates. However, we stress that the 
information needed for such a calculation is currently not up to date; therefore, 
specific quantitative results could change when constraints on the luminosity 
function are derived from the current pulsar sample. 

We calculate Pg(1Z) using the prior distributions of the two model pa- 
rameters for the pulsar luminosity function: the cut-off luminosity L m ; n and 
power-index p. We calculate these priors by fitting the marginal PDFs of 
L m i n and p presented by Cordes & Chernoff (1997). We obtain the follow- 
ing analytic formulae for /(L min ) and g(p): f(L min ) = a + aiL min + a 2 L 2 min 

and g(p) = lO /3o+/3lP+/32P , where a\ and 0\ (i = 0,1,2) are coefficients we ob- 
tain from the least-square fits and the functions are defined over the intervals 
£min = [0.0, 1.7] mJy kpc 2 and p = [1.4, 2.6]. We note that, although Cordes 
and Chernoff (1997) obtained /(L m in) over L m \ D ~ [0.3, 2] mJy kpc 2 centered at 
1.1 mJy kpc 2 , we consider f(L m \ n ) with a peak at ~ 0.8 mJy kpc 2 considering 
the discoveries of faint pulsars with L1400 below 1 mJy kpc 2 (Camilo 2003). 
We use the above priors to calculate P S (TZ): 



In Figure 2, we show the distributions of L m i n and p adopted (top panels) and 
the resulting global distribution of Galactic DNS merger rate estimates (bottom 
panel). We find that Pg(7Z) is strongly peaked at only around 15Myr — 1 . We 
note that this is a factor ~ 5.5 smaller than the revised rate from the reference 
model (TZ =83 Myr -1 ). At the 95% confidence interval, we find that the Galactic 
DNS merger rates lie in the range ~ 1-170 Myr -1 . These imply LIGO event rates 
in the range ~ (0.4 - 70) x lO^yr -1 (initial) and ~ 2 - SSOyr" 1 (advanced). 
Given these implications, it is clear that up-to-date constraints on L m \ n and p 
and their PDFs (a follow-up on Cordes & Chernoff 1997) are urgently needed. 

4. Rate constraints from Type Ib/c supernovae and binary evolution 



Based on our current understanding of DNS formation, the progenitor of the sec- 
ond neutron star is expected to form during a Type Ib/c supernova. Therefore, 
the empirical estimates for the Type Ib/c SN rate in our Galaxy can be used 
to provide upper limits on the DNS merger rate estimates. From Cappellaro, 
Evans, & Turatto (1999) we adopt ft SN ib/c - H00 ± 500 Myr^ 1 (for Sbc-Sd 
galaxies). Here, we assume H$ = 71km/s/Mpc and LB,gai = 9 x 10 9 Lb,suii- 

In order to find the fraction of SN Ib/c actually involved in the formation 
of DNS, we have examined population synthesis models calculated with the 
code StarTrack (Belczynski, Kalogera, & Bulik 2002; Belczynski et al. 2004) 
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Figure 2. The global P g (lZ) on a linear scale (lower panel) and the 
assumed intrinsic distributions for L m j n and p (upper panels). Dotted 
lines represent the lower (SN^) and upper (SN\j) bounds on the ob- 
served SN Ib/c rate scaled by 1/10 and 1/100 (see text). The empirical 
SN Ib/c rates range over ~ 600 — 1600 Myr -1 , where the average is at 
~4100 Myr -1 (Cappellaro, Evans, & Turatto 1999), beyond the range 
shown here. 
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and find very low rate ratios: 7 = ^/^-sNib/c ~ 0.001 — 0.005. Several models 
with He-star winds consistent with current observations (weaker than previously 
thought) lead to 7 ~ 0.005. We note that systematic overestimation of 7£gNib/c 
relative to SN II rates has already been pointed out (Belczynski, Kalogera, & 
Bulik 2002; this is related to the assumption of complete removal of H-rich 
envelopes). However, we find that this discrepancy would raise the value of 7 
by just a factor of a few. As an approximate constraint, we adopt the empirical 
^•SNib/c an d scale it by 1/10 and 1/100, reflecting the results from population 
synthesis calculations. We overlay these scaled values in Fig. 2 (dotted lines in 
the bottom panel) using the ranges for SN Ib/c reported by Cappellaro, Evans, 
k Turatto (1999). 

We note that our most optimistic DNS merger rate is 1Z = 224^^ Myr -1 
at a 95% confidence interval (Model 15 in KKL). We obtain 7 for SN Type 
Ib/c to be ~0.8 with the upper limit of 1Z at the 95% confidence interval. This 
corresponds to 7 ~0.1 with a SN Type II rate, which is factor 6.1 larger than 
that of SN Type Ib/c. In both cases, the most optimistic model is lower than 
the current empirical supernova rate estimates, but not really consistent with 
the results of population synthesis calculations. If we consider the global distri- 
bution, with the upper limit of 1Z at the 95% confidence interval, we obtain 7 
~0.15 and 0.025 for SN Type Ib/c and II, respectively. 



5. Prediction for more DNS detections by the PMB survey 

Acceleration searches of the PMB-survey data (Faulkner et al. 2003) should 
significantly improve the detection efficiency of DNS binaries. Although the 
data analysis is on-going, acceleration searches already led to the discovery of 
PSR J1756— 2251 (see Faulkner et al. 2004 and the contribition by Lyne in this 
volume) . 

Following the method described in Kalogera, Kim & Lorimer (2003), we 
calculate the probability to detect a pulsar similar to any of the observed DNS 
systems. We assume that acceleration searches can perfectly correct for the 
reduction in flux due to Doppler smearing, namely no degradation in the cal- 
culation of signal-to-noise ratio for the PMB survey is included. Considering 
observed DNS systems individually, we calculate the expected number of pul- 
sars to be detected by the PMB survey (N exp ). 

The probability distribution of A^ xp for each DNS pulsar sub-population i 
(B1913+16, B1534+12, J0737-3039) is given by: 



(3j 2 (iV exp + 1) 
(1 + Pi) 2 (1 + ' 



w-») = r, l :\ ,^ :z:^L > (2) 



where the constants are a measure of how less likely it is to detect pulsars 
without acceleration searches relative to with acceleration searches. For each 
sub-population, we calculate the mean values of N exp , which are iV exPi i9i3 = 0.9, 
^Vcx P ,i534 = 1-2, and iY cxPi0 737 = 1-9. 

The increase of the observed sample is very important for the reduction 
of the uncertainties associated with the DNS merger rate estimates. We note, 
however, that the discovery of new systems that are similar to the three already 
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Figure 3. P{N cxp ) for close DNS systems in the PMB survey. We 
consider all three observed systems in the Galactic disk. The mean 
value is -/V exp =4.0. The result shown here is based on our reference 
model. 

known does not necessarily imply a significant increase in the rate estimates. 
Significant changes are expected when new systems are discovered with pulse 
profiles or binary properties significantly different than the old ones, as those 
systems will reveal a new DNS sub-population in the Galaxy. 
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